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We present models where 7 + and 7 -, the exponents of the susceptibility in the high and low 
temperature phases, are generically different. In these models, continuous symmetries are explicitly 
broken down by discrete anisotropies that are irrelevant in the renormalization-group sense. The 
Zq-invariant models are the simplest examples for two-component order parameters [N = 2) and 
the icosahedral symmetry for N = 3. We compute accurately 7 + — 7 _ as well as the ratio v/v' of 
the exponents of the two correlation lengths present for T < T^. 
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The question of the equality of the critical exponents 
on the two sides of a second order phase transition has, 
apparently, not been raised for decades. The general 
renormalization group (RG) argument “showing” their 
equality goes as follows: Correlation functions are regu¬ 
lar in the presence of an external field, which is sufficient 
to go continuously from one phase to the other. More¬ 
over, if these functions satisfy the same RG equations 
above and below the critical temperature Tc, the same 
should hold true for the scaling behavior of quantities 
such as the susceptibility, the correlation length or the 
specific heat. Since the renormalization properties of a 
theory are identical in its symmetric and spontaneously 
broken phases, it follows that the critical exponents are 
identical in both phases (see for instance [1-3]). This is 
indeed what happens generically. 

To the best of our knowledge, Nelson [4] was the first 
to propose a counterexample based on the 0(2) model 
in dimension d = 3, to which is added either a cubic 
(CA) [4-12] or hexagonal (HA) anisotropy [4-6]. These 
anisotropies are taken into account in the Ginzburg- 
Landau hamiltonian by terms of order 4 and 6, respec¬ 
tively, which are irrelevant in the RG sense at the tran¬ 
sition. The corresponding fixed point is thus 0(2) sym¬ 
metric. However, Nelson argued [4] that they are danger¬ 
ously irrelevant [6, 13] in the low-temperature phase and 
that they, therefore, induce a modification of the expo¬ 
nent 7_ of the susceptibility. A rather counterintuitive 
result is that the difference 7+ — 7_ is larger for HA than 
for GA, whereas HA is “more irrelevant” than GA. A de¬ 
tailed study of the literature shows that, up until now, 
this striking result has been completely ignored. 

Because of its relationship with either deconfined quan¬ 
tum critical points [14] or pyrochlore [15] and the possible 
existence of two distinct phase transitions [16], the three- 
dimensional XY model with HA (and more generally the 
Zq-invariant models) has been studied again [17, 18]. Al¬ 
though only one transition has been found, the Zg models 
were shown to exhibit two correlation lengths below Tc, ^ 
and , that scale with two different critical exponents, iz 


and iz'. All authors agree that v/v' depends on the scal¬ 
ing dimension of the irrelevant HA term, but there are no 
less than three different scaling relations predicting this 
ratio, as well as several values obtained by Monte Carlo 
simulations [14, 17, 18]. 

In this Letter, we present a mechanism, valid not only 
for the XY case, to generate different critical exponents 
above and below Tc- The mechanism relies on the possi¬ 
bility of explicitly breaking a continuous symmetry down 
to a discrete one by terms that are irrelevant in the RG 
sense. In addition to the Z,j-invariant models, we build an 
example for Heisenberg spins showing icosahedral sym¬ 
metry. Using the nonperturbative renormalization group 
(NPRG), it is—contrary to perturbation theory—easy to 
show that the exponents 7-1- and 7_ are generically differ¬ 
ent for these models and easy to compute them, as well as 
v/v', accurately. Our approach allows us to completely 
clarify the physics of these models. 

Let us discuss the general idea underlying the differ¬ 
ence between 7+ and 7_. For concreteness, we con¬ 
sider a XY or Heisenberg model described by an 0{N)- 
invariant Hamiltonian (iLo(v))i to which is added a dis¬ 
crete anisotropy term t{x), H = i7o(Af)+'''an '’'(x)- We 

assume that t{x) is irrelevant in d = 3. The fixed point 
(FP) describing the phase transition is therefore 0{N) 
symmetric = 0)- If this term were irrelevant in the 
ordinary sense—that is, could be neglected (Aan = 0)— 
the model would be identical to the 0{N) model. It is 
important to remember that, in this case, not only the 
transverse (xt) but also the longitudinal (xl) suscepti¬ 
bilities diverge for all T < T), because of the Goldstone 
modes [19-21]. However, since the symmetry is discrete 
when Aan 7^ 0, there are no Goldstone modes and the 
susceptibilities cannot diverge for T < Tc- Thus, Xt^ 
vanishes only at Tc, and its scaling with AT = T — Tc 
obviously depends on the way Aan scales to zero close 
to the fixed point. Since this scaling is given by the 
scaling dimension of t{x), the exponent yx defined by 
Xx^ ~ {Tc — T^"^ for T < Tc cannot be equal to 7+. 
This is why t{x) is said to be dangerously irrelevant for 
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T < Tc ■ The same holds true for Xl^- 

Let us now give two examples, with N = 2 and 3, of the 
kind of anisotropy that produces a difference between 7 _|_ 
and ^T,L- We choose discrete subgroups G of either 0(2) 
or 0(3). These subgroups must satisfy two constraints. 
First, in order to have only one phase transition, there 
must exist only one invariant quadratic polynomial of 
G. Therefore, it must be z = 1, • • • , iV as in the 

0{N) model. Second, the interaction term that explicitly 
breaks the 0{N) symmetry must be irrelevant. A term of 
order 4 can be irrelevant compared to the 0(fV)-invariant 
term {ipi(pi)‘^. For N = 2 and d = 3, this is the case, 
for instance, of the term (pf + of cubic anisotropy. 
For reasons that are explained below, this kind of terms, 
being “weakly irrelevant”, induces only small differences 
between 7 _|_ and Jt,l that, moreover, require very large 
systems to be observable. We therefore prefer to consider 
terms that are “strongly irrelevant”, because they are of 
degrees higher than 4. 

For iV = 2, all the Z^-invariant models (g-state clock 
models) with q > A satisfy the two conditions above 
because the first invariant polynomial in ipi that is not 
0(2)-symmetric is of degree q. For instance, for Zg, this 
invariant reads : t = + 6 (pi(p 2 ~ 20 (pf(p 2 [6] [22]. 

For N = 3, the situation is more constrained because 
only the icosahedral group satisfies the two conditions 
above [23]. For all the other discrete subgroups of 0(3), 
the first invariant polynomials that are not 0(3) symmet¬ 
ric are of degree 4 and, therefore, are at best only weakly 
irrelevant. For the icosahedral symmetry, the first non- 
0(3)-symmetric invariant polynomial is of degree 6 and 
reads [23, 24]: 

T =(4$ - 2){p'l - ip2){^2 - V’DiV’l - ‘pD + 22{pip2P3f 

+ (Pl + + pi)(pi + P2~^ 

( 1 ) 

where 4) is the golden ratio [25]. 

The NPRG is based on Wilson’s idea of integrating 
fluctuations step by step [26]. In its modern version, it 
is implemented on the Gibbs free energy F [27-29]. A 
one-parameter family of models indexed by a scale k is 
thus dehned such that only the rapid fluctuations, with 
wavenumbers jgl > fc, are summed over in the partition 
function The decoupling of the slow modes (|g| < k) 
in Zk is performed by adding to the original hamiltonian 
H a quadratic (“mass-like”) term which is nonvanishing 
only for these modes: 

Zk[J] = j D(pe^^{-H[ip\-/\Hu[p\+J-p) ( 2 ) 

with AiLfc]^)] = \ J^Rk{q^)pi{q)pi{-q)—vf\ieve, for in¬ 
stance, Rk{q^) = Zk^k"^ — q^)9{k‘^ — q^), with 9 the step 
function and Z^. the field renormalization constant—and 
J ■ p = Ji(x)pi(x). The coarse-grained Gibbs free en¬ 
ergy rk[<j>] is defined as the (slightly modified) Legendre 


transform of log Zfc[Ji]: 

rfc[<^]+logZfe[J] =J-(j) - ^ J^Rkiq‘^)M9)^ii-<l) (3) 

where (/)i(a;) is the thermal average of Pi{x). When k is of 
the order of the inverse lattice spacing A, all fluctuations 
in Zk are frozen by the Rk term and the mean-field ap¬ 
proximation becomes exact. With the definition (3), this 
implies that Ffe=A[<^] = H[4i] [27]. Since Rk=o{q^) = 0, 
rfc=o[<^] = r[^] and is, thus, the free energy that we want 
to compute. 

The exact flow equation of F^ reads [27] (see Supple¬ 
mental Material for more details): 

dt^k[4>] = ^^4dtRk{q'^){^k\q,-q-,4>] + Rkiq))~^] (4) 

where t = log(fc/A), Tr stands for an integral over q 
and a trace over group indices, and F), ’[q,—q-,<j)] is the 
matrix of the Fourier transforms of the second functional 
derivatives of rfc[(^] with respect to (j)i{x) and 

Since it is impossible to solve Eq. (4) exactly, we must 
make use of approximations. To capture the critical 
physics, the simplest nonperturbative approximation is 
the derivative expansion [30, 31]. We use the (improved) 
lowest order, the local potential approximation prime 
(LPA’), which consists of retaining only a potential term 
in rfc[(^] together with a field renormalization constant 
Zk in front of the kinetic term [30, 31], 

rr^'['A] = ^ |^Z,[V0,(a:)]2 + UkiM^))^ ■ (5) 

The running potential Uk is defined by ^Uk{4>i) = 
T k [4’i ]) where the fields (jji are constant and is the vol¬ 
ume of the system. Its flow is obtained from Eq. (4) 
where F), Ms computed from (5) and is then evaluated in 
a constant field configuration (pi. 

On top of the LPA’, Ukipi) can be expanded around 
one of its running minima 0™“, which corresponds, at 
fc = 0, to the stable state of the system when Jj = 0, 

Uk = ^—(p-Kk) +Uoi^kT+ ^^{p- Kk) H- (6) 

with p = (j)i4>il2 and Kk = </>™</>™/2. The flows of 
Umn,k and Kk are obtained from that of Uk by acting 
with dt on both sides of their definition, 

dUk, , ,,, 

Umn,k - QpmQ.j-n \p=x.k,r=0 , \p=Kk,T=0 — U. J 

We have performed the calculations up to order 12 in 
the field-expansion equation (6). However, for the sake 
of simplicity, we present below the RG flow obtained for 
the XY model with HA within the simplest ansatz that 
includes only U 2 o,k and itoi.fc, which we call u and Xq 
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(we omit the k index in the following to alleviate the no¬ 
tation). In this case, corresponds to one of the six 
minima of Uk ; we call the direction pointing towards 
longitudinal and the perpendicular direction transverse. 
Once diagonalized, splits as usual into the two in¬ 
verse longitudinal and transverse propagators that each 
depend on a (running) “mass” 2uk and ISAe^^, which we 
call mn and tot- Finding RG fixed points requires us to 
work with dimensionless and renormalized variables. We 
thus define k = Zkk‘^~'^K, u = KdZ^^k’^~^u, and 

Ag = KjZ-;:^k^<^-^Xe, where = dT{d/2)2‘^-^Tr^ has 
been included for convenience. The flow equations read 
(see Supplemental Material for more details): 


dtk =(2 - d - rik)k + 


(1 ISitAg^ 

(2 + ^) 


hirnr) 



+ \h{rni) ( 8 a) 

dtu ={d — 4 + 2r]k)u — ISXehi'm^) + 9 {(^/ 3 (to|) 

-|-(m-I- 36KA6)^/3(mf-) (8b) 


dtXe ={2d- 6 + 3r]k)Xe 
-I-15A6(w + bicAg) 


l 2 {m\) - himl) 

"2 "2 
mi — 


( 8 c) 


FIG. 1. XY model with HA in d = 3. (a) Flow for T < Tc 
of the running anomalous dimension, (b) Projections of RG 
trajectories in the u — X% plane. G, XY and NG are the Gaus¬ 
sian, critical 0(2) and low-temperature (Nambu-Goldstone) 
fixed points, (c) Flows for T < Tc of the inverse transverse 
and longitudinal susceptibilities. The dashed line represents 
the flow of in the pure 0(2) model, (d) Flow of u for 
T < Tc showing its five different regimes. The two plateaus 
correspond to the vicinity of the XY and NG FP. 


with Inix) = 2(1 -I- a;)“”(l — ^/(d-b 2)), to| = 2uk, 
rn^ = 18A6K^. The running anomalous dimension is de¬ 
fined by 77fc = —dt log Zk which tends, at criticality, to 
the anomalous dimension ry for fc —)■ 0 [30, 31]. We show 
its flow in Fig. 1(a). 

The flow equations (8) are very simple, though they 
are nonperturbative. They show two crucial features. 
First, they take automatically into account the role of the 
masses tol,t and their decoupling: as long as fc > niL,T) 
that is, 1 > tol,T) the contributions coming from /„(tol) 
(resp. /„(to^)) are non vanishing; they become negligible 
when fc ^ tol,T) and the longitudinal (resp. transverse) 
mode is then said to decouple from the flow. Second, 
once generalized to arbitrary N, they reproduce the low- 
T expansion of the 0{N) nonlinear-sigma model at one 
loop in d = 2 -b e (when Ag = 0) and are also one-loop 
exact in d = 4 —e [31, 32]. Gathering all these properties 
in a single set of flow equations is out of reach of the 
usual perturbative expansions. 

Equations (8) admit three fixed points (with noninfi¬ 
nite couplings u and Ag) as shown in Fig. 1. The NG FP 
has coordinates = 2 — d/2, Ag = 0, = 00 . By 

integrating Eqs. (8) in d = 3 with different initial con¬ 
ditions, we find the RG trajectories shown in Eig. 1(b). 
Eor T < Tc, Ag and k reach fixed, nonvanishing values 
for fc < This is why the dimenionful inverse trans¬ 
verse susceptibility Xt^ stops running beyond this scale, 
see Eig. 1(c). The dimensionless analogues of these cou¬ 
plings, Ag and k, keep running according to their canon¬ 


ical dimension. Thus, for fc ^ 
k{k) ^ fc(r') - «XY 

(9) 

Ag(fc) ~ Ag(r') (fcc)""-® ~ Airr'"«' 

where ye is the scaling exponent of Ag and A|f is the 
initial value of Ag, that is, its value at the scale of the 
inverse lattice spacing. As for u, we find in d = 3, for 
T slightly below T^, five different regimes represented on 
Eig. 1(d). In region (v), u diverges and u reaches a finite 
value. 

The values of fc at which the RG flow departs, respec¬ 
tively, from the XY and NG FP define two length scales 
called ^ and 5', see Fig. 1. The first one, is the Joseph- 
son length (the correlation length of the amplitude mode) 
of the pure 0(2) model [33] [the anisotropy plays no role 
in part (ii) of the flow, see Fig. 1(d), when ^ is large be¬ 
cause Ag(fc ^ ^“^) ^ 1]. As long as the RG trajectory 
remains close to NG, the flows of u and k remain similar 
to the flows of the 0(2) model in the low-T phase [region 
(iv) in Fig. 1(d)], and xj} decreases as it does in the 
pure 0(2) model where the fluctuations of the Goldstone 
modes make it vanish, see Fig. 1(c). However, because 
Ag 7 ^ 0, the flow departs from the NG FP [14, 17]; u{k) 
stops flowing, as does see Fig. 1(c). This occurs 
when m^{k = — 1, which defines . This correla¬ 

tion length diverges at Tc according to ^ AT~‘' , by 
definition of v'. Using the definition of to|- and Eqs. (9), 
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it is now straightforward to show that 

u'= u{l + \y^\/2). ( 10 ) 

This relation was already obtained in [18], although we 
do not fully agree with its derivation [34]. The deriva¬ 
tion above shows that v' is universal, contrary to what is 
suggested in [15]. 

The scaling relations among exponents are readily de¬ 
rived from the discussion above and the usual scaling 
behavior of the potential, 

U = s^U n(s), Xe{s), ■■■ , (11) 

where s is a rescaling factor, the dots stand for the infinite 
set of irrelevant couplings, and Z{s) is the field renor¬ 
malization factor. Z{s) ^ g-d+ 2 -rixY^ ^ (^A)“^, 

where yyxv is the anomalous dimension at the XY FP, 
and Z{s) ~ for s <C (^A)“^ since rjk vanishes away 

from it (in particular at the NG FP), see Fig. 1(a). In 
the low-T phase, by taking two derivatives of (11) with 
respect to (j) in either in the transverse or longitudinal 
directions, and then taking (f> at the minimum of U, we 
obtain 

oc s‘^Z{s)Xe{s)k^{s) ( 12 ) 

oc s‘^Z(s)m(s)k(s). (13) 

By taking s ^ (C^)“^ in Eq. (12) we obtain: 'yx = 'y+ + 
I'lyel, and by taking s ~ (^'A)“^ and using Eqs. (9), (10) 
we obtain xl = 7-i- + (4 — d)v\ye\/2 [35]. Notice that 
the scaling relations derived above for Zg are generically 
valid. 

We have computed yq up to order 12 in field-expansion 
equation ( 6 ) to obtain converged results, see Table I. We 
observe that our value of ly' for Zg is very close to the one 
deduced from the scaling law Eq. (10) and Monte Carlo 
simulations in [18]. This validates our approach. We find, 
of course, that \yq\ increases with g; we thus deduce the 
rather counterintuitive result that the more irrelevant the 
anisotropy term, the larger the difference between 7 + and 
7 t.l- However, since diverges extremely rapidly close 
to T~ for large \yq\, it must be difficult to observe the 
scaling behavior of xl in a finite-size system for “large” 
values of q. As for xt, its measurement should not be 
more difficult than in the pure 0(2) model. Reciprocally, 
if \yq\ is too small, the transient regime before reaching 
the XY FP is very large [region (i) in Fig. 1(d)], and, 
thus, the corrections to scaling are also large; this spoils 
an accurate determination of the leading scaling behavior 
in finite-size systems. This is probably the case of CA in 
d = 3, where we find |j/ 4 | = 0.042 [at six loops |j/ 4 | = 
0.103(8) [12]]. 

We have presented a general mechanism leading to a 
large and measurable difference between critical expo¬ 
nents in the high- and low-T phases and a theoretical 
approach to compute them. For the XY case, we have 


Symmetry 

^4 

Zs 

Ze 

Zio 

y' 

0.71 

1.06 

1.44 

3.84 


0.72 [14] 

1.05 [14] 

1.6 [14] 
1.45 [18] 

2.8 [14] 

+ 

1 

0.029 

0.06 [ 12 ] 

0.74 

1.49 

1.58 [4] 

6.29 


TABLE I. Critical exponents in d = 3 for the XY, models. 
In both cases, the first row corresponds to our results. For 
the icosahedral symmetry {N = 3), we find at order eight in 
field expansion Eq. ( 6 ): i/' = 1.51 and 7 t — 7 + = 1.54 . 


resolved the existing discrepancies between the results 
obtained in Z^-invariant models [14, 15, 17, 18]. Let 
us also emphasize that layered decagonal quasicrystals 
[36, 37] showing tenfold anisotropies and XY spin sys¬ 
tems with HA [4] exist, which would enable a direct mea¬ 
surement of 7t,l ~ 7 - 1 - and v'. Another very interesting 
challenge is the possibility of measuring susceptibilities in 
Heisenberg systems with icosahedral anisotropy, possibly 
in quasicrystals. We recall that, for TV = 3, there are 
problaby many other anisotropies that are dangerously 
irrelevant; this likely would lead to differences between 
'yT,L and 7 + that are smaller than in the icosadrehal 
case, but that are possibly also measurable. Finally, it 
would be extremely interesting to investigate the two- 
dimensional [38-47] case with the NPRG approach. At 
the price of avoiding any field truncation and working 
in at least the second order of the derivative expansion 
[48-50], this is reachable. We leave it for future work. 

We acknowledge D.R. Nelson for pointing out to us 
the existence of quasicrystals with tenfold anisotropies, 
P. Sindzingre for discussions and clarifications about nu¬ 
merical issues, A. Rangon for pointing out to us some im¬ 
portant references, and F. Benitez, H. Ghate, N. Dupuis, 
F. Rose, M. Tissier and N. Wschebor for discussions and 
suggestions regarding the manuscript. 
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